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Abstract 

Spontaneous nucleation and the consequent penetration of vortices into thin 
superconducting films and wires, subjected to a magnetic field, can be consid- 
ered as a nonlinear stage of primary instability of the current-carrying super- 
conducting state. The development of the instability leads to the formation of 
a chain of vortices in strips and helicoidal vortex lines in wires. The boundary 
of instability was obtained analytically. The nonlinear stage was investigated 

by simulations of the time-dependent generalized Ginzburg-Landau equation. 
74.60.Ge, 74.76.-w 



1 



Typeset using REVT^ 



Finite resistance in massive type-II superconductors subjected to a supercritical magnetic 
field Hci is caused by the motion of vortices The order parameter is suppressed 

at a distance whose magnitude comparable to the correlation length ^ within the core of 
the vortex, and the induced currents and fields decay at the penetration length A. The 
dynamics of vortices both in massive three-dimensional systems and in two-dimensional 
films has became the subject of intensive studies (both analytical p and numerical [§]). 

Another limiting case is quasi-one-dimensional systems (transverse dimensions smaller 
than ^ and A) such as narrow microstrips and channels. Their resistance is caused by phase- 
slippage centers (PSC), which nucleate spontaneously in some narrow region of currents 
even without a magnetic field 

We consider here systems with transverse geometry intermediate between the two cases 
described above; i.e., for which the radius R of 3D wires satisfies the condition: ^ <^ R <^ \. 
Likewise for strips one assumes that hd ^ A^ together with h ^ C,, where h is the thickness 
and d is the width of the strip. We assume that the homogeneous external magnetic field is 
applied normal to the strip and parallel to the wire's axis. The current is applied lengthwise 
to the strip and parallel to the wire's axis. Our aim is to elucidate the onset and penetration 
of vortices into a sample. 

The problems of the surface instabilities and consequent nucleation of the vortices in 
bulk superconductors was widely discussed (see, e.g. Ijl^l)- The critical magnetic field for 
bulk semiinfinite superconducting sample, known as the " superheating field" , was obtained 
by several authors (see, e.g. [Jill). However, to our knowledge, for thin superconducting 



films and wires with intermediate geometry {hd <^ A^ or ^ <^ i? ^ A) such problems were 
not studied. 

One can expect that the stationary, spatially inhomogeneous current-carrying supercon- 
ducting state becomes unstable when the applied current j and/or magnetic field B exceed 
some critical value(s). Because of the translation invariance along the longitudial direction, 
the perturbations arising for the order parameter have the form of a standing wave with 
finite modulation wavenumber q. In virtue of the inhomogeneity of the stationary state, the 



perturbative eigenfunctions should be mostly localized near the edges. Exponential growth 
in time of the perturbations finally results in the formation of single zeroes (lines of zeroes 
for wires) of the order parameters on a set of points placed at the edges of a strip at a 
distance 27r/g apart. Such zeroes serve as the "seeds" for vortices which will be driven away 
later on from the edges by the current and fields. In this article we present the analytic 
treatment of this problem. Our results are supported by the numerical simulations of the 
TDGLE in the region where the analytic approach fails. 

Although the vortices in thin strips and wires are not topologically different from the 
usual vortices in massive superconductors, the intermediate geometry results in significant 
differences. In contrast to massive superconductors the magnetic field remains practically 
homogeneous on the scale of the strip or wire. As a result, only the vortex phase rather 
than the magnetic flux carried by the vortex is quantized Moreover, the values of critical 
currents and fields obtained happens to be strongly dependent on the transverse dimension. 
The vortices in thin strips and wires are different also from the PSC nucleating in a very 
thin wire (i? ~ ^ ) where one needs the potential difference across a sample, and instability 
typically appears only in a narrow band of currents 0. Note that still another type of 
vortices occurs in superfluid jets above the critical velocity |]T2[, which formally corresponds 
to zero magnetic fleld with the superconducting current greater than the depairing current. 

We start from the celebrated time- dependent Ginzburg - Landau Equation (TDGLE) in 
the form 0] 

u{dt + i/i)^ = (V - 2iAf ^ + (l - 1^1^) ^, (1) 
j = (Vv9 - 2A) - (V/i + 2dtA) , (2) 
V-j = ,V-A = 0, (3) 

where is the (complex) order parameter, cp = arg"^, A and fi are vector and scalar 
potentials, and j is the current density. The value of the parameter u has to be obtained 
from microscopic theory [Q. The usual dimensionless units are used 0. 

Eqs. dll-^) apply for thin strips and wires. Indeed, the condition dh <^ for strips and 
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i? ^ A for wires enable us to neglect the magnetic field created by currents ||I3| , [T6| and omit 
thereby the appropriate Maxwell equation. 

Let us start with strips. We choose the origin of the coordinate frame at the mid-point 
of the strip with the x axis lengthwise and the y axis in the lateral direction, so that the 
edges are located at {x,—d/2) and {x,d/2). The normal to the strip magnetic field B is 
described by the vector potential A = {—By, 0,0) (see Fig.l ). 

The TDGLE has a stationary current-carrying superconducting solution of the form 

y) = F{y) exp{ikx), /i = (4) 
jiy) = FiyY{k + 2By) (5) 

where F{y) satisfies the following nonlinear equation 

dlF+ [l-F^ -{k + 2ByY)F = (6) 

subject to the no-flux (superconductor-vacuum) boundary conditions dyF{d/2) = 
dyF{—d/2) = 0. This solution is the only stable superconducting state for B < Be. 

For the case of weak magnetic fields {B <^ k/d) one can find a perturbative solution of 
Eq. (I). In lowest order the solutions for F{y) and the mean current jo = 2 I~d/2j(y)^y 
have the form 

p^p^_^J^^V^kBs.nHV2F,y) 

Fo F2 cosh(V2Fod/2) ^ ' ^ ^ 



where jo = A;(l - k"^) - kd^{B^ + ©(l/d^)) and Fo = VT^^ (note that F{y) - Fq is 
antisymmetric only for Bd ^1). For arbitrary values of the magnetic field B Eq. (|^) can 
be solved only numerically. For the numerical solution of the nonlinear boundary problem 



we used the deferred-difference method [|I4|. The characteristic form of F(y) is shown in 
Fig. 2. 

We seek the perturbative solution of the 

TDGLE in the form \1/ = (^F{y) + ({x,y,t)^ exp{ikx). Linearizing TDGLE with respect 
to (, /i, sphtting real and imaginary parts of ( = d + ib, and representing the solution in 



the form (a, b, jj) ~ (a.(y), b{y), ^{y)) exp{iqx + ut) , we obtain the following equations for the 
functions a{y), b{y), fi{y) 

cuua = d^a + {Z{y, q) - 2F^)a - 2iq{k + 2By)b 
ojub = dyb + Z{y, q)b + 2iq{k + 2By)a - uFfi 
Fuuib = d^fi - {q^ + uF^)fi (8) 

where Z{y, q) = 1 — {k + 2ByY — F^ — q^. For convenience (in order to deal only with real 
coefficients in the Eqs. we redefine b ib. Then the Eqs. (P) are of the form 

uua = d^a + {Z{y, q) - 2F^)a - 2q{k + 2By)b 
ujub = dlb + Z{y, q)b - 2q{k + 2By)a - uF^i 
Fujub = d^fx - (g2 + uF^)fx (9) 

Eqs. represent the eigenvalue problem for u which has to be solved for the no-fiux 
boundary conditions for a, b and /i. The stationary solution becomes unstable when for the 
first time the negative uj{q) achieves zero. For u = the equation for fi has no solutions 
obeying the boundary condition. Therefore, we are left with two equations for a and b 
(^ = 0) 

d^a + {Z{y, q) - 2F^)a - 2q{k + 2By)b = 

dlb + Z{y,q)b-2q{k + 2By)a = (10) 

The system (|TD|) in general requires numerical solution. These equations contain 
four control parameters, namely, the magnetic field B, the mean current jo which is 
determined by k, the width of the strip d and the modulation wavenumber q. The 
numerical solution was implemented in the following way: 1) The stationary solution 
F{y) was determined (numerically) and this was then substituted into Eqs. (0); 2) 
For fixed q, Eqs. ( p!0D were solved for two sets of boundary conditions at the lower 
edge: dyai{-d/2) = dybi{~d/2) = 0,ai(-rf/2) = l,6i(-rf/2) = and dya2i-d/2) = 
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dyb2{—d/2) = 0,a2{—d/2) = 0,b2{—d/2) = 1; 3) Then the solutions at the second edge 
were examined; 4) The eigenvalue q was obtained from the condition 

dyaidyb2 - dya2dybi = 0. (11) 

The critical curves for the onset of instability in the (j, B) plane for several widths of the 
strip d are depicted on Fig. 3. The eigenfunctions a{y),b{y) are shown in Fig. 2. It can be 
seen that for B ^ the eigenfunctions are strongly localized near the edge of the strip where 
the suppression of the order parameter is maximal. The critical curves are terminated at the 
depairing current jp = 2/^/27 ^ 0.387, which corresponds to kp = l/VS. In this region a 
very weak magnetic field is needed to destroy superconductivity. In the practically relevant 
limit d ^ 1 (well-separated edges), Eqs. (|1^) can be essentially simplified by adiabatic 
elimination of a{y) from the second Eqs. (p!OD resulting in: 

^;'--^^a- ,_,t^%'|;)r,,. )^+w. (12) 

For k kp also g ^ and the critical curve is given by the solvability condition J_[i/2{^ — 
2-2{k+2By)^ )'^^ ~ ^' "^^^ Straightforward calculations results in the following conditions: 
jp — j = kjjB^d'^ and B^d'^ = —36/kp where 6 = k — kp <^ 1. However this trivial solution is 
relevant only for very small 6 ~ l/d"^. For larger 6 this curve merges with another curve which 
turns out to be unstable for smaller B. This curve obeys the scaling {q'^d, 6d, Bd"^) ~ 0(1). 
For this case the Eq.(0) is reduced to the Airy equation: 

d'^b=ia + /3y)b (13) 

where a = ^g^(— 16-\/35(l — g^) + — 725^ — 3g^), = —6\^q'^B (these expressions were 
obtained with the aid of the Maple program for analytical calculations). The critical curve 
for the Eq. (O) is given by the conditions (cf. ([Tl|)) 



dyAl{-fl)dyBl{-f2) - dyBl{j,)dyAl{j2) = 
^ [dyAt{^,)dyBt{^2) - dyBt{^,)dyAt{^2)] = (M) 



where 71^2 = /?^'^^(tc^/2 + a/ (3). The critical curves obtained by solution of Eq. (|T^) are 
shown also on Fig. 3. The critical curves obtained from the Airy equation (|13|,|1^ coincide 
with those within the thickness of the line. We conclude that the "d ^ 1" approximation 
works extremely well. As can be seen from Fig. 3, the critical field Be decreases with d, and, 
finally B^. vanishes for d ^ 00. This is consistent with well-known results [Q that the critical 
magnetic field in an infinite film is always zero. Figures 2-3 support the physical picture of 
an appearance of a chain of vortices at one edge of the strip (where the eigenfunctions are 
localized). It is useful to mention that the critical curves appear to be insensitive to the 
particular properties of the superconductor given by the material constant u. 

In order to follow the consequent stages of the instability, which cannot be covered 
by the linear analysis, we performed direct numerical simulations of TDGLE (with u = 
5.79 [^). The initial conditions were chosen typically as a homogeneous superconducting 
state perturbed by small amphtude noise. We used no-flux boundary conditions in the 
transverse direction (9y\E' = 0) and normal metal- superconductor boundary conditions in 
the longitudial direction (\l/(x,?/) for x — 0, L, where L is the strip length). The 



numerical scheme was the generalization of the split-step method described in |Ta,0,17 



the number of the grid points was 256x 128 and the timestep was 0.05 — 0.1. Results of the 
simulations, shown in Fig. 3, 4 clearly agree with the stability analysis. 

The analytical solution and the numerical simulations of TDGLE clarify all stages of 
the process: initial development of small perturbations due to applied current and/or field 
(Fig. 4a), their localization at some points at the lower edge (Fig. 4b) and an appearance 
of the chain of single vortices and its consequent tearing off and propagation inside the strip 
(Fig. 4c). Then the process repeats itself. When the vortex is sufficiently far away from 



the edges, its motion can be described, in principle, by a simple first order equation |15,16 



This periodic nucleation of vortices at the edge was observed for sufficiently wide range of 
transport currents < j < jp. Thus, one has a dynamic resistive state near the critical 
line. A gradual increas of the field B leads to faster nucleation, also the vortex motion 
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itself looses the regularity and reminds of some sort of vortex turbulence [1T5|,|6|. In that case 



one has also nonperiodic oscillations of the voltage across the sample, as it was observed in 
numerical simulations pr5|,|6|. Very close to the critical current jp the film doesn't exhibit 



dynamic resistive states and undergoes direct transition to normal state. Alternatively, for 
j ^ the period of the nucleation diverges. 

Our analysis can be generalized for the case of superconducting wires when B is parallel 
to the axis z of wire. Then the potential can be chosen in the form A = {Ar, Ae, A^) = 
(0,Sr/2,0). 

The equation for the stationary order parameter has the form (compare with Eq. (^)) 
d^F + ^drF +{l-F^-e- (Brf) F = (15) 

Eq. (|15D is subjected to the boundary conditions drF{0) = drF{R) = 0. 

We seek a perturbative solution in the form \E' = (-F(r) + ({r,6, z,t)) exp{ikz) Notice, 
that unlike Eqs. @, one now has the full 3D problem. In contrast to strips the wires 
also possess rotation symmetry. Therefore, the perturbative solution should break both 
translation and rotation symmetry, creating, thereby, helicoidal structure. 

Linearizing now the TDGLE, and substituting the perturbations in the form (a, b, fi) ~ 
exp(m^^ + iqz + ut) (a, b have the same meaning as for a strip), we are left for the critical 
value uj = with the following equations 

A„a„ + (Zi(r, q) - 2F^)an - 2qkbn + 2Bnbn = 

Anbn + Zi{r, q)bn - 2qkan + 2Bnan = (16) 



where Zi(r, q) = I - P - [Erf - F^ - q\ A„ = + \dr - ^. 

The boundary conditions read dr-an = drbn = for r = R and a„, 6„ ~ r" for r — > 0, 7^ 
0. For n = the conditions translate to dra^lO) = 9r&o(0) = 0. The numerical solution of 
Eqs. (|T6]) is similar to that of the strip if the fact that they contain the rotation number 
n is taken into account. Critical curves for each n intersect in some complicated way on 
the (j, B) plane. In order to select the most unstable rotation mode we used the following 
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algorithm. We reproduced the critical curves for the 6 first modes (n = 0..5, other modes 
turned out to be unimportant for the chosen radius of the wire) as a function of k, and 
then for each given k we selected the mode which is unstable for lowest value of 5. As 
a result, the critical curve in the (j, B) plane consists of several pieces of different n. As 
can be seen from Fig. 5 the rotation symmetry changes with B, namely, for small B the 
axi-symmetric mode n = appears while for larger B the higher helocoid solutions occur. 
The number of the "important" modes grows with the radius of the wire. For R ^ 1 the 
critical curves can be described by the reduction to effective Airy equation as that for the 



strip. For macroscopic wires the helicoid structures are well-known (see, e.g., [Q). 

It is a challenging problem to generalize the above formalism to a magnetic field orthog- 
onal to the axis of wire. However, the stationary state is then described by a 2D nonlinear 
partial differential equation, and the analysis is very hard. One can speculate that the 
critical curves for this case will not much differ from that of a strip with d ~ 2R. 

Finally, we have described all the stages of the appearance and the propagation of a chain 
of vortices in thin strips and wires, in contrast to the usual approach where the existence 
of vortices is assumed ad hoc. All existing methods of experimental observation of fluxon 
dynamics, e.g., voltage-current characteristics, temporal pulses at a measured voltage, etc. 
may be used to check our analysis. High-temperature superconductors are most suitable for 
these experiments. Indeed, they have large values of A and, hence, samples can be prepared 
satisfying the conditions dh ^ and R <^ X. In particular, for the YBaCuO system with 
parameters A = 1500 A, c/ = 500 A and /i = 10 A, one has dh/X^ <^ 1. For this system the 
scale of physical units shown in Fig. 2-5 are: i/phys ~ 1/ x lOA, Bphys ~ i? x 5 ■ lO'^ Gs and 
jphys = J X 2.5 ■ 10® A/cm^. The characteristic time scale is given by tphys = t x 10~^ sec. 
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FIGURES 



FIG. 1. Coordinate frame for the superconducting strip. 

FIG. 2. The amphtude of order parameter F{y) and the perturbative solutions a{y),b{y) 
as functions of y obtained from numerical solutions of Eq. (^) and Eqs. (^) for q = 0.2804, 
d = 30,k = 0.3, B = 0.014 and j = 0.22. 

FIG. 3. The critical Qc and Be versus j for two different values of d. The diamonds and heavy 
dots show the loci of stable and unstable stationary solutions of TDGLE obtained by numerical 
simulations. Dashed lines show the critical curves are given by Eq. (|l^). 

FIG. 4. Results of numerical simulations of the TDGLE. The time runs from the top to the 
bottom: a) t = 90, h) t = 150 and c) t = 210. The parameters used are d = 30,B = 0.018, j = 0.2 
and the length of the strip L = 70. The current is applied along x-axis and the magnetic field is 
perpendicular to the strip. On the gray coded snapshots of |^'(x,y)| ( = is shown in 

black and = 1 is shown in white) one clearly sees the onset and propagation of vortices. 

FIG. 5. The (j, B) critical curve for a wire with R = 10 
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